(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 9.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       157,          7]
NotebookDataLength[     12254,        394]
NotebookOptionsPosition[     11235,        354]
NotebookOutlinePosition[     11573,        369]
CellTagsIndexPosition[     11530,        366]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[BoxData[
 RowBox[{"\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"eta0", "=", 
    RowBox[{
     RowBox[{"1199169832", "/", "10000000"}], "Pi"}]}], "\[IndentingNewLine]", 
   RowBox[{"Dout", "=", "1"}]}]}]], "Input",
 CellChangeTimes->{{3.567937831575729*^9, 3.567937836394609*^9}, {
  3.567938028325918*^9, 3.567938053952236*^9}, {3.5689665795019817`*^9, 
  3.568966621318551*^9}}],

Cell[BoxData[
 FractionBox[
  RowBox[{"149896229", " ", "\[Pi]"}], "1250000"]], "Output",
 CellChangeTimes->{3.567938223207395*^9, 3.568966643175906*^9}],

Cell[BoxData["1"], "Output",
 CellChangeTimes->{3.567938223207395*^9, 3.5689666431778097`*^9}]
}, Closed]],

Cell[BoxData[
 RowBox[{
  RowBox[{"Z0", "[", "m_", "]"}], ":=", 
  RowBox[{
   RowBox[{"eta0", "/", "8"}], "  ", 
   RowBox[{
    RowBox[{"EllipticK", "[", 
     RowBox[{"1", "-", "m"}], "]"}], "/", 
    RowBox[{"EllipticK", "[", "m", "]"}]}]}]}]], "Input",
 CellChangeTimes->{{3.567937807009571*^9, 3.567937822628722*^9}, {
  3.567938066584167*^9, 3.5679381973386393`*^9}, {3.5689759707355757`*^9, 
  3.568976021106731*^9}}],

Cell[BoxData[{
 RowBox[{
  RowBox[{"LR", "[", 
   RowBox[{"m_", ",", "ratio_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"EllipticK", "[", 
     RowBox[{"1", "-", "m"}], "]"}], "/", 
    RowBox[{"EllipticK", "[", "m", "]"}]}], "-", 
   "ratio"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"M", "[", "m_", "]"}], ":=", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"Sqrt", "[", 
       RowBox[{"1", "-", "m"}], "]"}], "-", 
      RowBox[{"Sqrt", "[", "m", "]"}]}], ")"}], "^", "4"}], "/", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"Sqrt", "[", 
       RowBox[{"1", "-", "m"}], "]"}], "+", 
      RowBox[{"Sqrt", "[", "m", "]"}]}], ")"}], "^", "4"}]}]}]}], "Input",
 CellChangeTimes->{{3.568976073073997*^9, 3.568976108475854*^9}, {
  3.56897631536998*^9, 3.568976372675528*^9}}],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"sol", "=", 
  RowBox[{"FindRoot", "[", 
   RowBox[{
    RowBox[{"LR", "[", 
     RowBox[{"x", ",", 
      RowBox[{"16", "/", "17"}]}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"x", ",", "0.00001", ",", ".99999"}], "}"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"N", "[", 
  RowBox[{
   RowBox[{"Z0", "[", 
    RowBox[{"M", "[", 
     RowBox[{"x", "/.", "sol"}], "]"}], "]"}], ",", "16"}], "]"}]}], "Input",
 CellChangeTimes->{{3.5689761257217007`*^9, 3.568976175705579*^9}, {
  3.568976207749597*^9, 3.5689762566582117`*^9}, {3.568976383170842*^9, 
  3.5689763936214952`*^9}, {3.568976428878067*^9, 3.568976447821085*^9}, {
  3.5689767209946423`*^9, 3.568976723795498*^9}, {3.568976758282834*^9, 
  3.568976761580203*^9}, {3.568977476105802*^9, 3.568977536218994*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"x", "\[Rule]", "0.5660258102832174`"}], "}"}]], "Output",
 CellChangeTimes->{{3.568976219898591*^9, 3.568976258449953*^9}, {
   3.568976394559165*^9, 3.568976448644556*^9}, 3.568976726299491*^9, {
   3.568976761990006*^9, 3.568976773555294*^9}, {3.568977477003592*^9, 
   3.568977536577897*^9}}],

Cell[BoxData["204.24657010915863`"], "Output",
 CellChangeTimes->{{3.568976219898591*^9, 3.568976258449953*^9}, {
   3.568976394559165*^9, 3.568976448644556*^9}, 3.568976726299491*^9, {
   3.568976761990006*^9, 3.568976773555294*^9}, {3.568977477003592*^9, 
   3.5689775365806093`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"data", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"{", " ", 
     RowBox[{
      RowBox[{"di", "/", "Dout"}], " ", ",", 
      RowBox[{"N", "[", 
       RowBox[{
        RowBox[{"Z0", "[", 
         RowBox[{"di", "/", "Dout"}], "]"}], ",", "16"}], "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"di", ",", 
      RowBox[{"1", "/", "10"}], ",", 
      RowBox[{"9", "/", "10"}], ",", 
      RowBox[{"1", "/", "10"}]}], "}"}]}], " ", "]"}]}]], "Input",
 CellChangeTimes->CompressedData["
1:eJxTTMoPSmViYGAQAWIQvfHux2Lfqa8cE024y0H0nMOr6kH0zoJfPSD68Am/
ySC6+e6POSD6klHrIhBd1PEVTMcx+C0H0Txz9qwG0d9Mz3gWAOkgfycvEF01
vSgQRDeskIoD0ZN01ZJA9IbllskgOrhdLB1EK534C6YZrjgWgmiJrtllIFpD
0nk6iLa45b4ERD/yNt4Mol09osC0wdnFu0H09wmee0D0iv/yzc7TXznmGVV2
gGgO7Wk9IJrvZedsEH2B/8B8EH2jzmshiJZyaNgIopOy+zaB6Ft2/adAdNuV
PadBtMtz95sguirGB0w3JWr5egDpRwrN4SAaALoqpTY=
  "]],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     FractionBox["1", "10"], ",", 
     "75.2930835876484027458850068824808995564`16."}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["1", "5"], ",", 
     "64.04747972919402577918744441150425279323`16."}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["3", "10"], ",", 
     "57.0232377981369280571784587454309705356`16."}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["2", "5"], ",", 
     "51.64931542231276019283507305690297163254`16."}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["1", "2"], ",", 
     "47.09128918272133193352480005253991353359`16."}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["3", "5"], ",", 
     "42.9355064778395391990466722506810882642`16."}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["7", "10"], ",", 
     "38.88922485848638630253470891007135836412`16."}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["4", "5"], ",", 
     "34.62414955697107337904383442462680278888`16."}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["9", "10"], ",", 
     "29.45276526374695889116123048428645872215`16."}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.5680101021394444`*^9, {3.5680101394030657`*^9, 3.568010165777265*^9}, 
   3.5680104313749733`*^9, 3.568010461459507*^9, {3.568010493352314*^9, 
   3.5680105313899717`*^9}, 3.568010651893381*^9, 3.568010724545289*^9, 
   3.568966877639402*^9, {3.568966911015706*^9, 3.56896692308739*^9}, {
   3.568967054240767*^9, 3.568967061276733*^9}, 3.568967263960381*^9, {
   3.568967359898954*^9, 3.568967371458042*^9}, 3.5689765734856853`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"data", "=", " ", 
  RowBox[{"Prepend", "[", 
   RowBox[{"data", ",", 
    RowBox[{"{", 
     RowBox[{"\"\<d/D\>\"", ",", "\"\<Sqrt(e_r)Z0\>\""}], "}"}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.568966951261189*^9, 3.56896698827803*^9}, {
  3.568976585196211*^9, 3.5689765867227*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"\<\"d/D\"\>", ",", "\<\"Sqrt(e_r)Z0\"\>"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["1", "10"], ",", 
     "75.2930835876484027458850068824808995564`16."}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["1", "5"], ",", 
     "64.04747972919402577918744441150425279323`16."}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["3", "10"], ",", 
     "57.0232377981369280571784587454309705356`16."}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["2", "5"], ",", 
     "51.64931542231276019283507305690297163254`16."}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["1", "2"], ",", 
     "47.09128918272133193352480005253991353359`16."}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["3", "5"], ",", 
     "42.9355064778395391990466722506810882642`16."}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["7", "10"], ",", 
     "38.88922485848638630253470891007135836412`16."}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["4", "5"], ",", 
     "34.62414955697107337904383442462680278888`16."}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     FractionBox["9", "10"], ",", 
     "29.45276526374695889116123048428645872215`16."}], "}"}]}], 
  "}"}]], "Output",
 CellChangeTimes->{
  3.568966995147019*^9, 3.568967161145849*^9, 3.568967270447404*^9, {
   3.568967363125724*^9, 3.568967374171385*^9}, 3.568976591221856*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Grid", "[", 
  RowBox[{"data", ",", 
   RowBox[{"Frame", "\[Rule]", "All"}]}], "]"}]], "Input"],

Cell[BoxData[
 TagBox[GridBox[{
    {"\<\"d/D\"\>", "\<\"Sqrt(e_r)Z0\"\>"},
    {
     FractionBox["1", "10"], "75.2930835876484027458850068824808995564`16."},
    {
     FractionBox["1", "5"], "64.04747972919402577918744441150425279323`16."},
    {
     FractionBox["3", "10"], "57.0232377981369280571784587454309705356`16."},
    {
     FractionBox["2", "5"], "51.64931542231276019283507305690297163254`16."},
    {
     FractionBox["1", "2"], "47.09128918272133193352480005253991353359`16."},
    {
     FractionBox["3", "5"], "42.9355064778395391990466722506810882642`16."},
    {
     FractionBox["7", "10"], 
     "38.88922485848638630253470891007135836412`16."},
    {
     FractionBox["4", "5"], "34.62414955697107337904383442462680278888`16."},
    {
     FractionBox["9", "10"], "29.45276526374695889116123048428645872215`16."}
   },
   AutoDelete->False,
   GridBoxFrame->{"Columns" -> {{True}}, "Rows" -> {{True}}},
   GridBoxItemSize->{"Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}],
  "Grid"]], "Output",
 CellChangeTimes->{3.568967008818071*^9, 3.568967163511722*^9, 
  3.568967273263843*^9, 3.568967376720477*^9, 3.56897659466486*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Expand", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"Sqrt", "[", 
       RowBox[{"1", "-", "m"}], "]"}], "-", 
      RowBox[{"Sqrt", "[", "m", "]"}]}], ")"}], "^", "4"}], "/", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"Sqrt", "[", 
       RowBox[{"1", "-", "m"}], "]"}], "+", 
      RowBox[{"Sqrt", "[", "m", "]"}]}], ")"}], "^", "4"}]}], "]"}]], "Input",\

 CellChangeTimes->{{3.5689768595155563`*^9, 3.568976890000218*^9}}],

Cell[BoxData[
 RowBox[{
  FractionBox["1", 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      SqrtBox[
       RowBox[{"1", "-", "m"}]], "+", 
      SqrtBox["m"]}], ")"}], "4"]], "-", 
  FractionBox[
   RowBox[{"4", " ", 
    SqrtBox[
     RowBox[{"1", "-", "m"}]], " ", 
    SqrtBox["m"]}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      SqrtBox[
       RowBox[{"1", "-", "m"}]], "+", 
      SqrtBox["m"]}], ")"}], "4"]], "+", 
  FractionBox[
   RowBox[{"4", " ", "m"}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      SqrtBox[
       RowBox[{"1", "-", "m"}]], "+", 
      SqrtBox["m"]}], ")"}], "4"]], "-", 
  FractionBox[
   RowBox[{"4", " ", 
    SuperscriptBox["m", "2"]}], 
   SuperscriptBox[
    RowBox[{"(", 
     RowBox[{
      SqrtBox[
       RowBox[{"1", "-", "m"}]], "+", 
      SqrtBox["m"]}], ")"}], "4"]]}]], "Output",
 CellChangeTimes->{{3.568976868452402*^9, 3.568976890907304*^9}}]
}, Open  ]]
},
WindowSize->{982, 1037},
WindowMargins->{{36, Automatic}, {43, Automatic}},
FrontEndVersion->"9.0 for Linux x86 (64-bit) (November 20, 2012)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[579, 22, 386, 9, 77, "Input"],
Cell[968, 33, 153, 3, 49, "Output"],
Cell[1124, 38, 94, 1, 32, "Output"]
}, Closed]],
Cell[1233, 42, 425, 11, 26, "Input"],
Cell[1661, 55, 831, 26, 55, "Input"],
Cell[CellGroupData[{
Cell[2517, 85, 804, 19, 55, "Input"],
Cell[3324, 106, 335, 6, 32, "Output"],
Cell[3662, 114, 286, 4, 65, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[3985, 123, 884, 23, 32, "Input"],
Cell[4872, 148, 1712, 46, 129, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6621, 199, 313, 8, 32, "Input"],
Cell[6937, 209, 1502, 44, 129, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8476, 258, 119, 3, 32, "Input"],
Cell[8598, 263, 1157, 28, 300, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9792, 296, 498, 16, 32, "Input"],
Cell[10293, 314, 926, 37, 69, "Output"]
}, Open  ]]
}
]
*)

(* End of internal cache information *)
